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THE CALCULATION OF ANTENNA RADIATION PATTERNS 
BY A VECTOR THEORY USING DIGITAL COMPUTERS 


Richard F. Schmidt 
Advanced Development Division 


ABSTRACT 


This document reviews the Kirchhoff -Huygens -Fresnel theory of diffrac- 
tion, including Kottler boundary correction, and discusses the reduction of this 
formulation for the special case of metallic reflectors with unbounded conduc- 
tivity. The formal transition from the scalar to the vector case is given in de- 
tail. Radial and transverse field components are identified, the entry of the 
near-field terms is discussed, and the physical quantities (sheet current, charge 
distribution, etc.) associated with the mathematical development are presented 
through a simple dimensional analysis. The analytic aspects of the problem are 
considered, and the points of departure from conventional analysis to numerical 
methods and digital computation are outlined. 
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GLOSSARY OF NOTATION 


Symbol 


Meaning 


vector operator del 


r, R 


radial distances, subscripted 


ip (x, y, z, t) scalar wave with space and time dependence 


normal to surface or line 


surface, subscripted 
volume, subscripted 


an infinitesimal value 


time 


frequency 
angular frequency 


x, y, z 


Cartesian coordinates 


wave number 


velocity of light 


electric and magnetic field vectors, subscripted 

electrical conductivity, subscripted 

magnetic permeability, subscripted 

inductive capacity, subscripted 

dielectric displacement and magnetic flux density 

surface charge density 


sheet current, subscripted 
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GLOSSARY OF NOTATION (Continued) 


Symbol 

C 

dl 

M, L, T, Q 


J 

A 

» # 

J 

0 




1 i, etc. 


i, j , k 


A 



Mean ing 

contour 

differential length, a vector 

fundamental quantities: mass, length, time, and volume change 
density 

current density 

delta, an increment indication 

a constant equal to j/-l in complex variable theory 
a scalar field, generic symbol 

corresponds to, goes over to, or approaches (as required by 
context) 

given that 

implies that 

dyads or second order tensors when vector dot or cross product 
operation is omitted 

unit vectors in the x, y, z Cartesian coordinate directions 
respectively 

a unit vector in the radial direction, spherical coordinate system 
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THE CALCULATION OF ANTENNA RADIATION PATTERNS 


BY A VECTOR THEORY USING DIGITAL COMPUTERS 


INTRODUCTION 

This document is the first in a scries of documents describing the develop- 
ment of a general purpose, modular computer program for the calculation of 
antenna radiation patterns. The program represents an outgrowth of several in- 
vestigations of monopulse antennas and the inertialess displacement of monopuise 
antenna patterns in space, carried out during the period 1966-1967, In the course 
of these investigations it was found that reliable volumetric radiation pattern in- 
formation was frequently unavailable and that direct measurement was both 
time-consuming and costly, Range and component imperfections added to the 
measurement problem. A transfer function was desired to bridge the gap be- 
tween prime -feed and secondary radiation patterns. 

Examination of the several approaches available for calculating antenna 
radiation patterns showed that the current distribution method (Ref, 1, p. 144} 
appeared to offer the best compromise between the extremes of complete rigor 
and the approximations of the simple geometrical optics approach. The current 
distribution method provides a vector solution to the problem. It leads to the 
satisfactory solution of real physical problems where the size of the scattering 
body is large compared to the operating wave length. The current distribution 
method supplies information on the amplitude, phase, and polarization of radia- 
tion patterns. It also yields radial as well as transverse field components and 

the electric and magnetic field integrals contain near-field terms proportional 
I 

to — r , permitting field computations in that region. Mutual coupling on the sent- 
ry 

tering surface and reaction on the prime radiator are not included. 

This document contains some material which is not readily available in pub- 
lished texts. Its main purpose is to review the vector theory principles required 
in the subsequent computations, including the current distribution method and the 
Kirchhoff-Huygens-Fresnel theory, using a common notation and convention 
system. It provides several frequently omitted steps in the development of the 
theory, and details the arguments which lead from a scalar to a vector formula- 
tion. The purpose of this first document is to provide a theoretical foundation 
for the analytical details and computational methods to be described in subsequent 
documents . 

Kirchhoff 's scalar theory is first reviewed and a vector analogue is written. 
The analogue is manipulated by means of vector identities until a convenient 
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form is obtained for calculating the scattered fields. Kottler’s correction term 
for open (finite) surfaces, is derived from boundary conditions, and several in- 
tegrals of the general formulation are set to zero after an examination of the 
boundary conditions established by the assumption of perfectly conducting sur- 
faces. Transverse and radial components of the final version of die scattered 
field equations are identified. The origin of induction terms (fields proportional 

1 

to — r) is discussed. Physical quantities such as sheet current, and charge dis- 

tribution are associated with each of the integrals comprising the general solu- 
tion. Finally, a brief discussion of analysis and numerical/computational 
methods is presented, 

Subsequent documents will present details relating to the source polariza- 
tion vector, field interaction at a boundary, feed rotation matrices, feed displace- 
ment vectors, surface normals, differential areas, and radii of curvature. Com- 
posite surfaces and distorted surfaces will also be considered, together with 
dual reflector systems. The function of various subroutines in a modular-type 
computer program will be discussed, and direct as weF n numerical integration 
techniques will be presented. Examples of radiation patterns obtained with IBM 
7094 and IBM 3G0 Mod 91 computers and SC 4020 plotters will be appended in 
these latter reports . 


KIRCHHOFF'S DIFFRACTION THEORY 

The integral theorem attributed to Kirchhoff is a straightforward mathe= 
matical description of the physical process described earlier by the Huygens - 
Fresnel principle. A detailed discussion can be found in the literature. The 
derivation of the scalar Kirchhoff integral theorem involves several assumptions 
and the limitations of the theory, which become evident in the derivation, are 
worth noting at the outset as these cause some difficulty later in the application 
of the vector formulation to real physical problems. 

A monochromatic scalar wave (\p) is assumed to originate from a region 
containing sources as shown in Figure 1. 

In a vacuum, the space dependent part of \p satisfies the time -independent 
scalar Helmholtz equation. 

(V 2 + k 2 )'Jj ~ 0 


where 


2 


V (x, v, z, t) ” v ('t/ y, *5 f “ ,a,t on surface S £ Sj + S 2 + S 3 



VOLUME (V) 

SURFACES (S V S 2 , S 3 , S') 


P (x‘, y', **) 


P (x, y, z) 


S 2 IS CLOSED AT INFINITY 


Figure 1, Diffracting System 


y2 


B 2 J2_ 

Br 2 r Br 


when \// = i/; (r) alone. 

k = o’/c = wave number 

o) - angular frequency of radiation 

c a velocity of light 

It is assumed that \f* has continuous first- and second-order partial derivatives 
inside and on surface S. If <// is any other function which satisfies the same con 
tinuity requirements, Gres Vs Theorem 2 holds and 

J v 2 - ^‘V 2 i/»)dv - ~ T^) ds = 0 

*Ref. 2, p. 375 and Ref. 3, p. 280. 

2 Ref. 4, p. 82, and Ref. 5, p* 22. 
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if (v 2 + k 2 ) 0' =0 also. Numerous functions satisfy the wave equation, but 
gikr e" ikr 

\p' s — ^ or — - — f conforms to the physics of the problem and satisfies 
the radiation condition, 

e ikr 

If [jj\ 5» « ■ ■■ — , then 4 ji not analytic due to the singularity at r = 0, and 

the continuity conditions for Green's Theorem are violated. The situation is 
corrected by surrounding the singularity with an infinitesimally small sphere S' 
of radius e for which n s r, Instead of | » 0, f + JJ, * 0 because of the 
singularity. 

The Kirchhoff Integral Theorem follows immediately. 


f 

s s' 





ds' 


s' as e 0 




e ik e 
6 


Bv// 

Be 


e 2 d 0 


H Solid AngJ e 

= + J" \fi d H = 4 tt \p (x 1 , y\ z' ) 

12 Solid Angle 
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It is noteworthy that the theorem was derived without any explicit reference 
to frequency or to the curvature of S, The description of the bounds t 'es S t and S 2 
was indefinite. In particular, the physical aspects of the boundaries vvere not 
discussed. It is possible to construct a vector analogue directly from the Kirch- 
hoff Integral Theorem above, however, it is instructive to follow the historical 
developments in optics to obtain a better insight into the physics of the problem 
before making this step. 

The Huygens -Fresnel principle (1818) included periodicity in both space and 
time, and represented the first significant explanation of diffraction phenomena. 
Heuristic arguments were invoked in the original derivations based upon this 
principle which was, in fact, detached from physical reality. Kirchhoff’s Integral 
Theorem (1882), as derived above, embodies the basic ideas of the Huygens- 
Fresnel principle and actually reduces to the latter for many special cases. The 
Kirchhoff Integral Theorem is superior since it is more inclusive and simpler. 
Many texts on optics now proceed with the derivation of the so called Fresnel- 
Kirchhoff diffraction formula, 


-i A [ [e ik ( r+s ) 

0 (K 1 * y ' > z ' ) ~ ~2XjJ ^ [ cos r ) ” cos (n, s )] ds , 

Where A is a constant, r is the distance from a source to an aperture point, s 
is the distance from aperture point to observer and ds is on the aperture alone . 
The physical picture is that of Figure 1. Good discussions on the limitations 
and assumptions of the diffraction formula can be found in the literature, and 
these are binding on the Kirchhoff Integral Theorem and its vector analogue. 
They are worth reviewing. 

In the derivation of the Kirchhoff diffraction formula, S 2 is usually taken to 
be a planar, opaque screen extending to infinity. If the problem is one in optics, 
Sj is large in terms of the operating wavelength. Kirchhoff originally assumed 
that 

3 1 b 

(1) 0—0 and = o on the interior of S 2 ; 

(2) the field at S t is identical with the incident field, and is unaffected by 
the presence of the opaque boundary. 

d\jj 

Fundamental analytic errors are involved in this procedure ! If \fj = = 0 

over S 2 , then there is a discontinuity about the contour c> which is the boundary 
between and S 2 , and Green’s Theorem is violated. The electromagnetic field 



cannot be represented by a single scalar wave function, but is characterized by 
a set of jsuch frmctions which gives the components of the electric and magnetic 
fields, E and H, Finally, Maxwell's equations should be satisfied, but the solu- 
tion is defective in this respect also, at this stage of the development, since_the 
field equation for E cannot be obtained from the curl of the expression for H , 
and conversely. The computation of \p behind S 2 , by Kirchhoff’s methods, can 
lead to non-zero values if the wavelength is sufficiently large. This is contrary 
to the initial assumption. In spite of these difficulties, the classical Kirchhoff 
theory leads to satisfactory solutions for many diffraction problems. 1 

The success of the Kirchhoff approach in physical optics problems is due to 
the fact that the size and radius of curvature of the aperture S 1 is large in terms 
of wavelengths, polarization is usually ignored, only power density is recorded 
(on film, usually), phase is ignored, and most of the diffracted radiation is thrown 

forward so that the assumption dj = = 0 is largely satisfied. The adaptation 

on 

of the theory to radio -frequencies must, therefore, be made with discretion and 
the original assumptions of the theory should be reviewed upon application to 
physical problems. 


VECTOR EXTENSION OF THE KIRCHHOFF INTEGRAL THEOREM 

The formal vector extension of the Kirchhoff Integral Theorem is obtained 
directly from the scalar theorem 


0(x', y', z') 




Let 


JL eikr 

477 r 


Then 


*Ref. 6, pp. 461-464 


6 



\p(x\ y 1 , z' ) 



{ 


[Gn • V0 - 0n • VG] dS 


which is still a scalar formulation. If E » 


E(x', y', z') 



C G(n • V)E - E(n * V)G] dS 


which is the vector formulation 1 for the electric field. E is a complex vector, 
representing the illumination of the aperture whose diffraction pattern is to be 
calculated. 


The vector transformations which lead to the field equations 2 are straight- 
forward, with the exception of a single dyad operation. A detailed derivation is 
included in Appendix A as the manipulations are somewhat tedious. The result 
is: 


E(x', y', z') - - [i coil (n x H)G + (n x E) xVG + (n • E) VG] dS , 

(closed) 

for the electric field and by similar process, or from the principle of duality 

H(x', y\ z') = + ^ [ioj e (n x E)G - (n x H.) x VG - (n • H) VG] dS 
(closed) 


for the magnetic field. 


DISCONTINUOUS SURFACE DISTRIBUTIONS 

The application of Green’s Theorem in the derivation of the scalar Kirchhoff 
Integral Theorem makes it necessary that the illumination i p be continuous and 
have continuous first- and second-order partial derivatives. All of the derivations 


^ef. 6, p. 464 and Ref. 3, p. 283 
2 Ref. 6, p. 469, equations (29) and (31) 
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up to this point of the development have tacitly assumed that this was so, and 
that the integration was over a closed surface^as in Figure X. The formal vector 
extension must now be modified since £ and H will be discontinuous at the edge 
of reflecting surfaces. Kottler has given a contour distribution consistent with 
the requirements of the problem. 1 A discontinuity in the tangential components 
of E and H in passing from a physical surface to a void implies an abrupt change 
in the surface current density. This can be accounted for, according to Kottler, 
by an accumulation of charge on the contour. 

Boundary conditions at the surface are 2 : 

(1) n. (B a - Bj) » 0 

(2) n . (D 2 - Dj) = 7) (reserving co for angular frequency) 

(3) n x (E 2 - Ej) = 0 

(4) nx (H 2 - H t ) = K. 

The subscripts identify the opposite sides of a boundary as shown in Figure 2. 



Figure 2. Contour at Edge of Reflecting Surface 


The figure shows region A as a void, and region S as a perfectly conducting 
surface , however , the dual of this could be assumed and would lead to the same 


^ef. 6, p. 468 

2 Ref. 6, p. 35 7] (surface charge density) and K (sheet current) can exist only when one cr(con 
ductivity) is unbounded. 
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basic arguments. Following' Kottler , 1 the discontinuity in the sheet current 
equals the line charge accumulation on the contour. 


n 2 * (K a - K s ) ~ i mV 

The K a and K s are obtainable via boundary relationships i 


K a s n x (H 2 - Hj) = n x (H. - H. ) = 0 , 


and 

K s = n x (H 2 - Hj) = n x (0 - 2H. ) = -2n x H. . 

Then 

•n 2 . (K a -K s ) = n 2 . (2n x H. ) = 2H. n 2 x n - 2H. -dlj = Hj • d fj , 

where 


n 2 x n 


dl 


l 


is in the direction of dl . 

Since the physical phenomena involves line charges, the radial vector VG is 

Hi • dl\ 

annexed directly to v ; , recalling the previous results for surface 

ICO 

1 

charge distributions. A factor — is also required for the development of an 

electric field in the present system of units. The total effect is obtained from a 
line integral along Contour C. The complete Kottler field correction terms are 
given in Stratton 2 as 


^ef. 7, p. 48 4 
2 Ref. 6, p. 468 
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E(x', y', Z') E - I jf ?G H, • dl 
and 

H(x', y', x<) ' r— * f VGE, • dT . 

V J ‘ ICUli J 0 1 

The latter integral is obtained by similar process or directly by the principle 
of duality. 

CALCULATION OF THE SCATTERED FIELD 

The results of the vector extension of Kirchhoff’s Integral Theorem are now 
recalled and joined to the Kottler contour integrals for discontinuous surface 
distributions. The surface integrals are taken over an open surface (Sj) and the 
hypothetical illuminating fields E and H are given subscripts (1) to emphasize 
the fact that the incident fields E i and H. have interacted with the metallic sur- 
face, Further notational changes are now introduced so that engineering forms 
will be available and comparison with available examples on applications will 
be facilitated. Using 


0 



the field equations satisfying Maxwell’s equations 


and 


V x E + 


BB 

Bt 


0 


VxH- 


BD 

Bt 


J 


become: 
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E(x' t y\ z') - 


1 JL_ 

j«?c 


J Vi/i H 1 


dl 


\V 

(Kottler) 


i 7 f 


4^ j ija^(n x Hi) ^ + (n x Ej) xVs// + (n ‘ E x ) V*/0 dS, 

8 i 


(Kirchhof f) 


H (x',y, !') 4^-/ W'Ei • cU 

c 

© 

(Kottier) 

+ 4 -" J Cj cue (n x E t ) i p - (ri xHJ xV 0 - (n * H x )V>p] dS , 

© © © 

(Kirchhoff) 


If the scattered field is to be calculated for metallic surfaces with unbounded 
conductivity, then boundary conditions can be examined for the surface integrals 
and some simplification results. From the preceding discussion the fields behind 
the physical surface vanish (E 2 = H 2 s 0). Also, e x = e 2 = 1, =_ /x 2 _= 1, _ 

'J\ = 0, o ~ 2 - co in the constitutive equations 1 D = e r e 0 E , B = fi r ii Q H, J = crE. 

Then, for this special case, the general equations can be simplified since 


n x (E 2 -Ej) = -n xEj - 0 


and 


n • 0 2 H 2 - - -/xpi * Hj = 0 

Surface integrals ® , (©> and ® can be set to zero. An examination of the line 
integral © shows that E t • dl refers to the tangential component of the total 


^ef. 6, p. 10 
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field in medium (1). The third boundary condition, n x (E 2 - Ej) = 0 states that 
the tangential electric fields are continuous at the boundary. jSince E 2 = 0, Uie 
tangential electric field in medium (1) must vanish. That is E t 0, but n x E t •*- 0 
and E t * dl = 0. By this argument integral © must also vanish. 1 

An examination of the operator (V) provides additional insight to the computa- 
tional aspects of the problem. Since 


i ft — 


p "* kr 

r 


is a function of (r) alone, take 


V«/f 


Bt// 


A 

Br 


Ml r 

Br 


2 


Then 


V\fj “ 




r 


for the near -field case and 


V0 =* - jk0 l r 

for the far-field case. Vector _Vy!/ is a local radial vector on dS. The radial 
and transverse components of E(x', y' , z') and H(x', y' , z') can be identified 
prior to calculation for the far-field case since the Cartesian components of all 
the I vectors are equal when the field point is taken at infinity, or at very large 
distances. 

For the distant E( x <, y',z') field, integrals Q) and (j) must be radial vectors 
due to the fact that the direction of Vt p is a unique’ radial vector. In general, in- 
tegral ® will contain transverse vector components (1 0 , 1^) and a radial (l r ) 
component. As the point of the observer (the point at which the pattern is 


^ef. 1, p. 147 

2 Ref. 4, p. 97; Ref. 8. Table II 
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calculated) moves toward Infinity the — terms become negligibly small 

r 2 

1 

with respect to the terms that diminish as . The summation of the radial 
components, proportional to for integrals (T)> (D> ancl ® can be expected to 
vanish also, leaving’ only transverse electric fields in the limit as r , 

For Hie H(x', y', z') field, only integral CD contributes to the solution, It is 
known that the cross-product of two vectors is a vector which is orthogonal to 
the given vectors. The direction of vector Vi p is a unique radial vector for the 
far-Eield case , and (n x H) x Vt/>is , therefore , orthogonal to I r , or transverse , 
Only transverse magnetic fields remain in die limit as r -» 

The general formulation, therefore, reduces to: 

EC*', y\ *') » ^~r~ | • <il - -±- / (n * H,)0 + (n • E,)V0] dS , 

(Kottler ) (Kirchhoff) 

H(x',y', z') = - J (n x H j)xV0 dS ; 

s 1 

(Kirchhoff) 


which is applicable to near- and far-field problems for perfectly conducting 
surfaces. 


ANALYSIS AND NUMERICAL METHODS 

The application of the field equations for E(x*, y', z') and H(x', y*, z«) to a 
variety of antenna configurations presents several problems. These will be dis- 
cussed in detail in succeeding documents, and only the general approach will be 
outlined here. If the reflectors which are to be considered can be described by 
a parametric or other suitable representation, then the normals to the surface (n) 
and the differential areas (dS) can frequently be derived in a straightforward 
manner. The illuminating fields at the surface (E i , H. ) are calculated next. A 
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polarization vector is needed Cor the source, and six degrees of freedom (three 
due to translation, three due to rotation) can be visualized for the placement and 
orientation of that source. Space divergence from source to surface (proportional 

to -~) is included. 

t* 

The incident fields are allowed to interact with the^surface according to the 
boundary conditions, local normal (n), etc, and E, and il, are established for the 
evaluation of the integrals for E(x', y', z‘). Conventional analysis (vector analy- 
sis, complex algebra, differential geometry) appears adequate up to tills point, 
but the evaluation of the integrals for arbitrary antenna configurations is directed 
toward numerical methods. 

The difficulties associated with the integration phase of the problem are 
already evident, 1 


SUMMARY 

This document provides a theoretical basis for the diffraction-pattern cal- 
culations performed by the Antenna Systems Branch, Advanced Development 
Division of the Goddard Space Flight Center. The vector Kirchhoff theory, with 
Kottlor correction for finite (open) surfaces, was presented without reference to 
particular geometric forms or antenna configurations, Subsequent reports will 
present detailed derivations of subroutines for iT, dS, etc. When using ellipsoidal, 
paraboloidal, spheroidal, hyperboloidal f conical and other surfaces. Composite 
and distorted surfaces will be considered, together with dual reflector systems. 
The relation of the various subroutines to a modular-type computer program will 
be discussed in these reports. Direct and numerical integration techniques will 
be presented. Examples of radiation patterns obtained with the IBM 7094 and 
IBM 360 Mod, 91 computers will be appended in the form of a series of plots 
from an SC -4020. 
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APPENDIX A 


VECTOR TRANSFORMATIONS FOR THE KIRCHHOFF INTEGRAL 


E(x\y\0 ~ [G(n • V)E- E(n • V)G]dS 
(n • V) GE - G (n • V) E + E (n * V) G 
[ ] s (n * V) GE - 2E (n - V) G 


Using 

Ax(BxC) e (A • C) B - (A • B)C, 
n x (E xVG) = (n • VG) E - (n • E) VG 
VG x (nxE) = (VG • E) n ~ (VG • n) E. 

Then 

C 1 = (n * V) GE - nx(ExVG) - (n • E) VG + VG x (nxE) - (VG • E) n . 

Using 

Vx (0 A) = 0Vx A + V0 x A , 

Vx (GE) = GVxE+VGxE. 


i 
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Using 


V * (<£A) ~ 0V • A + V0 • A , 

V • (GE) = GV • E + VG • E , 
and 

GV • E = 0 

in the absence of free charge. 

Then 

C ] = (n • V) GE + nx(V x(GE)) -n (V * (GE)) 

- m: (G(V x E)) - (n • E)VG + VGx(n x E) . 

Using the integral relationships (Ref. 6, Appendix II) 

(1) f (V • A)dv = f A • H dS 

j v J s 

(2) f (V x A) dv = f (n x A) dS 

Jy J s 

(3) f V<^> d v = <T 0 n dS 

Jv Js 

and making the correspondences 

(1) A - (VGE) , (2). A - Vx (GE) , (3) 0 - V • (GE) 

obtain Jackson’s (Ref. 3,p. 284) result 
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(1) f (VGE) • n dS 

(2) i" nx (V x (GE)) 

(3) f n (V • (GE)) dS = f V(V • (GE)) dv , 

J s J v 

Apparently Jackson uses (n • V)GE = n * (VGE) where the latter is a dyadic 
form (second-order tensor). Relationships in Morse and Feshbach (Ref. 11, p. 65) 
and Coburn (Ref. 12, p, 46, p. 113) show that this is admissible. In general, the 
distributive rule of vector multiplication holds so that 



3F « - 3F y . - 3F„ 


n • (VF) = i • (ii) -j- x n, + i ■ (ij) n, + i • (ik) — 


+ j ' 00 n y + j • (ii) 37 n y + i ■ (i k ) 37 n y 


a * /S 


BF 


BF 


A SS > 


BF. 


+ k ' ^ ki ) Tz n * + k * 17 n z + k * ( kk ) Bz n 


Bz 


and 


n • (VF) - (n • V) F . 


Then 

f [(n • V)GE + nx ( Vx(GE)) -n (V • (GE))] dS = 


{ [V2 (GE) 


+ V x V x (GE) -V(V • (GE))] dv 


0 


since 


V x V x V = V(V*V) -V 2 V . 
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The vector extension may be rearranged to read as follows: 


E(x\ y>, z>) = - £[nx(G(^ x E)) + (FT • E)VG + (n x E) X Vg] dS, 

(closed) 

e + i k r 

The texts of Stratton and Jackson use g = — - — and the time -varying term e“ xwt 

j k r 

consistently. Silver uses g = ^ 7 - and e +jwt * The present development follows 
the former so that 1 


V x E 


Bt 


3 He~ Ia)t , - 


k 2 = co^ jj . £ , c - (m 0 e 0 )" 1/2 • 


Then 


E(x', y', z' ) - - I* [ iaj/j, (n x H) G + (n x E) x VG + (n • E) Vg] dS 

Js 

(closed) 

which is part of Stratton’s equation (29), 2 


1 Re£. 6, p. 23 

2 Ref. 6, p. 469 
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APPENDIX B 


UNITS, DIMENSIONS, AND PHYSICAL QUANTITIES 

The M.K.S., system of units (Ref. 6, Appendix I) will be used. Units and 
dimensions for the field equations are tabulated below. 


Table 1 

Units and Dimensions 


Quantity 

Symbol 

Dimensions 

M.K.S. Unit 

Charge 

q 

Q 

Coulomb 

Current 

I 

T" 1 Q 

Ampere 

Charge density 

p 

l _3 q 

Coulomb/cubic meter 

Current density 

J 


Ampere /square meter 

Conductivity 

■ 

cr 

M _1 L" 3 T q 2 

Mho/Meter 

Electric potential 

0 

ML 2 r 2 Q” 1 

Volt 

Electric field intensity 

E 

MLT" 2 Q' 1 

Volt/meter 

Dielectric displacement 

D 

L -2 Q 

Coulomb/square meter 

Inductive capacity 

€ 

M' l L“ 3 T 2 Q 2 

Farad/meter 

Magnetic flux 

0 

ML 2 T _1 Q“ 1 

Weber 

Flux density 

B 

MT‘ 1 Q“ 1 

Weber /square meter 

Magnetomotive force 

m.m. f. 

T l Q 

Ampere -turn 

Magnetic field intensity 

H 

L" 1 T" 1 Q 

Ampere -turn/meter 

Permeability 

P 

MLQ" 2 

Henry/meter 


’•0 










The eight integrals of the complete formulation for E(x',y',z') and 
H(x’, y 1 , z 1 ) are examined to verify that the dimensions are correct throughout. 
Certain factors are then singled out and associated with physical quantities for 
the case when cr l - 0, o 2 - a>. 

Integral CD 

- • dT -*■ (T)(ML 3 T" 2 Q* 2 ) (L -2 ) (L- l r‘Q) (L) =»mlt- 2 q-> 

Jc 

/H l • dl\ 

where y — - — /=> (t) (L" 1 !" 1 Q) (L) => q => electric charge. 


Integral (T) 

-■3- J jcu/x (n X Hj) i/idS => (T* 1 ) (MLQ‘ 2 ) (L‘‘r l Q) (L* 1 ) (L 2 ) =»MLT- 2 Q _1 

s i 

where (nxHj) => lT 1 T " 1 Q »=> electric sheet current k (Ref. 6, p. 37). 

K is defined as 


JA1 


A1 - 0 

j _ co 


0*2 *-♦ 00 , 


Integral (T) 



where (n xEj) 


(n x Ej) x Vt//dS 
MLT " 2 Q" 1 => 


=* (MLT -2 Q" 1 ) (L" 2 ) (L 2 ) 


magnetic sheet current K m 


*=> MLT" 2 Q" 1 
(Ref. 1, p. 68). 
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Integral 0 


£ J <*•*. 


) V^clS => (MLT“ 2 Q~ l )(L“ 2 )(L 2 ) MLT“ 2 Q” 1 


where (FT * Ej) => MLT , ' 2 Q“ 1 , but (n • D) - e (n * E) L" 2 Q =^> electric surface- 

charge density r) (Ref. 1, p. 07), and <j 2 - oo . 


Integral (s 


; r -\ 

■siJ 


Jwi TZ l 1 ’ dl ** ( T )^ rlL “ 1( 5 2 H L “ 2 ) (M Lr2 Q~ l ) (L) => L” 1 T " 1 Q 
(E. * dl) 

where ** (T l ) (MLT"' 2 Q" 1 ) (L) «?> ML 2 T~ 1 Q“ 1 ^magnetic charge. 


Integral © 

J jweCnxE^i/zdS => (T _1 ) (M' l L _ 3 T 2 Q 2 ) {MLT' 2 <r‘)(L' 1 )(L 2 ) ^l’V'q 
77 S 1 

where (nxE t ) = ?: >MLT" 2 Q'' 1 => magnetic sheet current K m as before. 

Integral © 

- ~ J (n x H,) x V^dS = 4 . (V 1 T” 1 Q) (L~ 2 ) (L 2 ) =4> L- l T _1 Q 

Sj 

where (nxl^)^ L" 1 T'' 1 Q => electric sheet current K as before. 

Integral © 

- 47 / (n • Hj) V^dS => (L^T^Q) (L' 2 ) (L 2 ) =* 1/ 1 T" 1 Q 
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where (ivHj) => L " 1 r l Q, but (H ‘ K m (n • H x ) => (MLQ" 2 ) (L" l T" J Q) => 
M T” 1 Q ” 1 => magnetic flux density. 
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